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FOREWORD 

This  report  descirbes  a  portion  of  the  results  obtained  on  NASA  Grant 
NSG  3044.  This  work  was  done  under  subcontract  to  the  University  of 
Illinois,  Urbana,  with  Prof.  S.S.  Wang  as  the  Principal  Investigator.  The 
prime  grantee  was  the  Massachusetts  Institute  of  Technology,  with  Prof. 

F.J.  McGarry  as  the  Principal  Investigator  and  Dr.  J.F.  Mandell  as  a  major 
participant.  The  NASA  -  LeRC  Project  Manager  was  Dr.  C.C.  Chamis. 

Efforts  in  this  project  are  primarily  directed  towards  the  development 
of  finite  element  analyses  for  the  study  of  flaw  growth  and  fracture  of 
fiber  composites.  This  report  presents  analytical  and  experimental  re¬ 
sults  for  delamination  crack  growth  in  fatigue.  The  work  is  a  follow-up 
to  an  earlier  interim  report  NASA  CR-135248.  H.T.  Wang  was  a  research 
assistant  at  M.I.T.  when  the  experimental  work  was  performed. 
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Abstract 

This  paper  presents,  an  investigation  of  interlaminar  crack  growth 
behavior  in  fiber-reinforced  composites  subjected  to  fatigue  loading. 

In  the  experimental  phase  of, the  study,  interlaminar  crack  propagation 
rates  and  mechanisms  were  determined  for  the  cases  of  various  geometries, 
laminate  parameters  and  cyclic  stress  levels.  An  advanced  singular  hybrid- 
stress  finite  element  method  was  used  in  conjunction  with  the  experimental  re¬ 
sults  to  examine  tte  local  crack-tip  behavior  and  to  characterize  the  crack 
propagation  during  fatigue.  Results  elucidate  the  basic  nature  of  the 
cyclic  delamination  damage,  and  relate  the  interlaminar  crack  growth  rate  to 
the  range  of  mixed-mode  crack-tip  stress  intensity  factors.  The  study 
provides  fundamental  insight  into  the  problem,  reveals  several  important 
features  of  the  interlaminar  fatigue  failure,  and  should  be  of  practical 
importance  in  selection,  testing  and  design  of  composite  materials. 
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INTRODUCTION 

Damage  tolerance  has  been  of  great  concern  for  the  structural  reli¬ 
ability  and  integrity  of  fiber-reinforced  composites  in  their  advanced 
engineering  applications.  Interlaminar  cracking,  sometimes  also  called 
delamination,  is  considered  one  of  the  most  frequently  encountered  types 
of  damage  occurring  in  the  material  subjected  to  fatigue  loading.  It  is 
commonly  initiated  at  geometric  discontinuities,  processing  and  manufac¬ 
turing  defects,  machining  and  service  flaws,  and  takes  the  form  of 
separation  of  plies  under  mechanical ,  thermal  and/or  environmental  loading. 

The  presence  and  growth  of  an  interlaminar  crack  may  result  in  (1)  signifi¬ 
cant  reduction  of  stiffness,  (2)  exposure  of  the  interior  to  adverse  environ¬ 
mental  attack,  and  the  degradation  of  the  composite,  which  may  lead  to 
the  final  failure  of  the  structure.  Thus,  understanding  the  interlaminar 
crack  growth  behavior  under  cyclic  stresses  is  of  fundamental  importance 
in  the  application  of  composite  materials. 

Laminated  composites  generally  consist  of  a  number  of  plies  with 
different  fiber  orientations  separated  by  resin-rich  interlaminar  regions 
as  shown  in  Fig.  1.  The  interlaminar  region  usually  has  a  dimension  of 
one  twentieth  to  one  tenth  of  the  individual  ply  thickness,  depending  on 
manufacturing  and  processing  methods.  Since  the  strength  is  usually  low 
and  stresses  are  high  in  the  Interlaminar  region,  almost  all  kinds  of  defects 
such  as  voids,  inclusions,  broken  plies,  etc.,  unavoidably  tend  to  grow  in 
an  interlaminar  mode  at  very  low  nominal  stress  levels  and  to  propagate 
along  preferential  paths.  This  phenomenon  is  unique  to  composites  and  is 
particularly  significant  when  the  materials  are  subjected  to  cyclic  loading. 
The  development  of  a  through-the-thlckness  crack  from  a  defect  or  flaw  in  a 


3 


real  composite  may  be  difficult  in  many  cases.  This  is  in  contrast  to  the 
flaw  growth  behavior  in  monolithic  metals  in  which  defects  usually  develop 
into  a  through-the-thickness  crack  during  fatigue. 

The  importance  of  the  fatigue  delamination  has  long  been  recognized, 
but  research  progress  has  been  slow,  since  this  kind  of  real-life  defect  is 
characteristically  difficult  to  study  due  to  the  inherent  complexities  of 
the  crack  geometry,  the  material  anisotropy  and  discontinuity,  and  the 
crack-tip  singularity.  Because  of  the  matrix-dominated  and  localized  frac¬ 
ture  behavior  of  the  delamination,  it  would  be  unwise  to  make  any  simpli¬ 
fication  which  fails  to  consider  the  important  resin-rich  interlaminar  region. 
Numerous  studies  have  been  conducted  to  investigate  the  initiation  of  an 
interlaminar  crack  in  composites  caused  by  geometric  discontinuities  [1-3]*, 
processing  and  manufacturing  defects  [4,  5]  .  Relatively  little  work  con¬ 
cerning  the  interlaminar  crack  growth  behavior  has  been  reported  [6,  7], 
especially  that  associated  with  in-service  and  manufacturing  defects  under 
repeated  stresses.  This  paper  presents  an  experimental  and  analytical  inves¬ 
tigation  of  the  crack  growth  behavior  of  an  interlaminar  crack  initiated 
from  a  service-induced  surface  notch  in  a  composite  laminate  subjected  to 
cyclic  loading.  Unidirectional  composites,  which  have  the  simplest  laminate 
configuration,  are  examined  first  to  explore  the  fundamental  nature  of  the 
interlaminar  crack.  Interlaminar  fatigue  failure  behavior  of  more  compli¬ 
cated  angle-ply  composites  is  then  studied  and  compared  with  the  results 
obtained  from  the  basic  unidirectional  fiber  configuration. 


*Numbers  in  brackets  designate  References  at  end  of  paper. 
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EXPERIMENTAL  PROGRAM 

One  of  the  central  concepts  in  this  experimental  program  is  to  develop 
a  suitable  test  technique  in  which  stable  interlaminar  crack  growth  can 
be  obtained  under  cyclic  fatigue  without  being  affected  by  the  complex  three- 
dimensional  edge  effect  [8].  Experiments  were  conducted  to  identify  the  influ¬ 
ential  geometric  and  laminate  parameters,  to  examine  the  interlaminar  failure 

mechanisms,  and  to  characterize  rationally  the  fatigue  damage  behavior  in 
the  composites.  They  also  helped  to  establish  a  realistic  mechanical  model 
for  later  analytical  study  and  provided  experimental  data  to  validate  the 
results . 

Materials  and  Specimen  Preparation 

Scotchply  Type  1003  E-glass/epoxy  (manufactured  by  3M  Company,  Min¬ 
neapolis,  Minn.)  was  used  in  the  fabrication  of  all  test  specimens.  The 
material  was  received  in  a  unidirectional  prepreg  form.  Test  composites 
were  made  by  a  compressive  molding  process  following  manufacturer's  instruc-' 
tions .  Specimens  used  in  the  study  were  eight-ply  composite  laminates  of 
various  fiber  orientations  and  stacking  sequences.  An  interlaminar  crack 
initiator  was  Introduced  in  the  form  of  a  surface  notch  penetrating  one  or 
more  plies  as  desired  by  cutting  individual  plies  with  a  razor  blade,  prior 
to  the  laminate  assembly.  The  specimen  had  a  gauge  length  of  six  Inches  and 
a  width  of  one  inch,  and  the  surface  notch  extended  completely  across  the 
specimen  width  as  shown  in  Fig.  2. 

Experimental  Procedures 

All  interlaminar  crack  growth  tests  were  conducted  at  room  temperature 
in  a  laboratory  environment.  Fatigue  experiments  were  run  on  an  Instron 
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Model  1211  d5niamlc  cycler .  Specimens  were  cycled  under  a  controlled 
sinusoidal  load-time  condition  with  a  minimum  tensile  stress  of  1.0  ksi 
and  some  maximum  tension  at  a  frequency  of  approximately  5  cycles  per 
second  (Fig.  3) .  No  significant  heating  of  the  specimen  was  observed  at 
this  frequency  during  the  test.  A  traveling  optical  microscope  mounted 
on  the  testing  facilities  was  used  periodically  to  monitor  the  interlaminar 
crack  extension  in  the  specimen  under  transmitted  light  and  to  measure  ac¬ 
curately  the  delamination  crack  length  after  a  given  number  of  cycles . 

Fatigue  crack  growth  data  were  obtained  for  several  cases  in  which  inter¬ 
laminar  cracks  were  initiated  at  different  depths  under  the  surface. 
Experimental  results,  expressed  by  delamination  length  versus  number  of 
loading  cycles,  are  reported  in  the  next  section. 

Advantages  of  Present  Test  Method 

The  present  experimental  method  has  several  advantages  particularly 
suitable  for  studying  the  cyclic  interlaminar  crack  growth  behavior  in 
composites.  As  mentioned  previously,  the  initiation  of  delamination  from 
geometric  discontinuities  is  complicated  by  the  three-dimensional  edge  effect. 
Thus,  the  growth  of  the  interlaminar  crack  may  be  expected  to  be  influenced 
also  by  the  edges  of  the  laminate.  Better  appreciation  of  the  nature  of 
the  problem  can  be  obtained  by  investigating  the  delamination  growth  by 
itself  without  being  affected  by  geometric  boundaries.  This  purpose  may 
be  achieved  in  the  current  experiments  by  using  specimens  wide  enough  to 
minimize  the  localized  edge  effect  oh  the  propagating  . 

interlaminar  crack  front.  The  experiment  can  also  provide  a  controlled, 
stable  delamination  growth  at  a  desired  depth  and  enables  the  crack  length 
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to  be  measured  accurately  and  conveniently.  Moreover,  the  specimen 
geometry  is  considered  relatively  easy  to  model  in  the  analytic  study, 
which  will  be  discussed  in  the  next  section. 

ANALYTICAL  APPROACH 

The  current  analytical  approach,  using  an  advanced  singular  finite- 
element  method  of  the  hybrid -stress  model,  has  a  capability  to  overcome 
the  aforementioned  difficulties  and  can  provide  a  numerically  exact  solu¬ 
tion  for  the  problem.  The  method,  pioneered  by  Plan  et  al.  [9],  gives  a 
rapid  rate  of  solution  convergence  and  enables  an  accurate  description  of 
the  crack-tip  response  by  introducing  a  crack-tip  superelement  including 
singular  and  higher  order  terms  based  on  the  complex  variable  formulation 
of  Muskhelishvili  stress  functions.  It  is  particularly  suitable  for  studying 
crack  problems  in  composites,  since  the  complex  crack  geometry,  laminate 
properties  and  traction  boundary  conditions  of  crack  surfaces  can  be  exactly 
modeled.  Additional  features  of  the  method  include  the  flexibility  of 
formulation,  the  selection  of  element  characteristics,  and  the  expedient 
achievement  of  interelement  compatibility. 

Mechanical  Model  and  Assumptions 

Microscopic  observation  in  the  experimental  study  indicated  that  inter¬ 
laminar  crack  growth  is  clearly  a  matrix-dominated  progressive  failure 
mechanism  occurring  in  a  resin-rich  interlaminar  region.  Thus,  the  composite 
laminate  is  modeled  as  an  assembly  of  anisotropic  plies  bonded  together  by 
thin  resin  interlayers  subjected  to  an  in-plane  cyclic  stress.  Each  fiber 
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reinforced  lamina  of  thickness  has  its  material  axis  oriented  from 

the  loading  direction  by  an  angle  0^.  The  resin-rich  interlayer  is  assumed 
to  have  a  uniform  thickness,  t^,  and  perfectly  bonded  with  the  plies  every¬ 
where  except  the  region  where  a  crack  is  initiated  from  the  surface  notch 
tip.  The  crack  geometry  and  laminate  configuration  are  conveniently  ex¬ 
pressed  by  . )»  where  the  double  slash  and  the  under¬ 

line  represent  the  location  of  the  interlaminar  crack  and  the  penetration 
depth  of  the  surface  notch,  respectively.  The  interlaminar  crack  is  modeled, 
for  simplicity  without  losing  generality,  as  a  flaw  completely  embedded  in 
a  resin-rich  interlayer  as  shown  in  Fig.  4.  Studies  on  other  extreme  cases 
such  as  the  crack  located  at  a  vanishing  interlayer  or  at  a  ply/ interlayer 
interface,  giving  a  more  complicated  oscillating  three-dimensional  singular 
stresses,  are  reported  elsewhere  [10]. 

Method  of  Analysis 

General  procedure  of  formulating  the  plane  crack  problem  in  the  composite 
is  described  in  this  section.  Details  of  the  formulation  for  the  crack-tip 
superelement  and  surrounding  non-singular  elements  have  been  given  elsewhere 
[11,  12]  and  are  not  repeated  here.  Briefly,  the  formulation  of  the  analysis 
is  based  on  the  variational  principle  of  modified  complementary  energy. 
Following  MuskheHshvlll' s  formulation  [13],  the  stress  and  displacement 
fields  may  be  expressed  in  terms  of  two  stress  functions,  ((>(z)  and  ipCz)  , 
of  a  complex  variable  z.  A  mapping  function  w(5)  is  then  introduced  to 
transfer  the  singular  crack  domain  to  a  ^“Plane  by  z  =  w(5)  =  where 
-  it/2  ^  arg  ^  ^  ir/2.  The  functions  (|)(5)  and  are  analytical 

and  interrelated  by  the  traction  boundary  condition  along  crack  surfaces. 
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In  the  crack-tip  superelement  formulation,  is  assumed  to  have  the  form 

(1) 


n 


<P(0  =  I  b. 


and  function  may  be  obtained  as 


41(0  =  -  I[b.(-1)^  +jjh.]^^ 


(2) 


where 


or 


b  .  =  6 ,  +  i  6  ,  .  (non-symmetric  case) 
J  3  n+j 


b ,  =  (symmetric  case) 


(3) 

(4) 


and  3s  are  real  constants  to  be  determined.  Using  Eqs.  1  and  2,  the  stress 
and  displacement  fields  may  be  expressed  in  terms  of  the  stress  coefficients 
3s.  The  boundary  traction  T  can  be  calculated  by 


T.  =  a. .  V. 

1  ij  3 


(5) 


Interpolating  boundary  displacements  u  in  terms  of  nodal  displacements  g 
and  inserting  a  and  T  into  the  variational  principle,  the  crack-tip  super- 
element  stiffness  matrix  may  be  obtained. 

Surrounding  non-singular  element  stiffness  matrices  are  formulated  by 
a  conventional  hybrid-stress  finite  element  method  through  a  minimum  com¬ 
plementary  energy  principle  [14].  The  anisotropic  elastic  properties  of 
each  glass /epoxy  lamina  are  considered  in  the  analysis  by  introducing  an 
appropriate  compliance  matrix  in  the  element  stiffness  formulation.  Since 
boundary  displacement  functions  are  independently  assumed  for  the  crack 
element  and  non-singular  elements,  interelement  compatibility  can  be  ensured 
by  a  suitable  choice  of  interpolation  functions. 

The  assembled  governing  equations  for  the  whole  system  may  be  written  as 


Kq  =  Q 


(6) 
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where  K  is  the  global  stiffness  matrix,  and  Q,  the  consistent  loading  vector. 
After  solving  q  from  Eq.  6  by  an  appropriate  solution  scheme,  the  stress  and 
strain  fields  can  be  determined.  Stress  intensity  factors  K^.  and  cor¬ 

responding  to  an  applied  stress  level  and  an  interlaminar  crack  length  may 
be  found  from  the  relation 


=  >/27 

(7) 

and 

•'ll  ■  «n-U 

(8) 

The  associated  elastic  energy  release  rates  and  can  be  obtained  by 

a  standard  procedure. 

The  accuracy  and  convergence  assessments  of  the  solutions  are  compli¬ 
cated  by  several  unusual  features  of  the  problem  and  of  the  method  of  analysis 
due  to  the  inherently  singular  nature  of  the  interlaminar  crack.  A  study  of 
the  accuracy  and  convergence  of  the  analysis,  and  solution  stability  has  been 
conducted-  by  testing  cases  for  which  independent  solutions  are  available. 
Excellent  agreements  with  existing  closed-form  solutions  were  obtained. 

The  results  indicated  that  accuracy  within  approximately  one  percent  of 
the  converged  solutions  of  and  K^.^.  can  be  achieved  by  the  optimum  mesh 
arrangements  used  in  the  current  study.  Details  of  the  study  were  reported 
in  [12]  . 


For  the  cases  of  cyclic  loading,  the  ranges  of  interlaminar  crack-tip 
stress  intensity  factors,  AK^  and  may  be  determined  by 


AK^  =  (K^)  -  (K^)  . 

I  I  max  I  min 


(9) 


and 


AK  ^  =  (K^^)  -  (K^,)  . 

II  II  max  II  min 


(10) 


where  (K  )  and  (K  ) 

I  max  II  max 


are  opening  and  shearing  mode  stress  intensity 
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factors  corresponding  to  a  maximum  cyclic  stress  level  a  for  the  delami- 

max 

nation  crack  of  length  £ , ,  and  (K^)  .  and  (K^^)  .  are  the  values  associated 

d  I  min  II  min 

with  a  minimum  cyclic  stress  level  a  .  .  Numerical  results  of  (K  ) 

min  I^max’ 

,  (K  )  .  and  (K  )  .  were  accurately  computed  for  the  delamination 
II  max  I  min  II  min 

test  specimens  by  the  computer-aided  numerical  method. 

RESULTS  AND  DISCUSSION 

Typical  results  of  the  crack  growth  study  of  the  eight-ply  fiberglass 
composites  subjected  to  nominal  cyclic  stresses  are  presented  in  this 
section.  In  the  experimental  phase,  interlaminar  crack  growth  rates  were 
determined  by  the  constant  load-amplitude  fatigue  tests.  In  the  analytical 
phase  of  the  study,  the  crack-tip  stress  field  and  the  ranges  of  the  cyclic 
stress  intensity  factors  were  obtained  for  the  laminates  with  various  geometric 
and  ply  '  configurations .  Interlaminar  crack  growth  rates  during  fatigue 
were  then  characterized  by  the  ranges  of  the  mixed-mode  cyclic  stress  intensity 
factors.  The  Individual  glass/epoxy  lamina  used  in  the  analytical  modeling 
had  a  thickness  dimension  of  0.01  inch,  and  the  interlaminar  layer  was  taken 
as  one-twentieth  of  the  ply  thickness  as  observed  in  the  microscope.  Elastic 
constants  for  each  lamina  were  given  by 

=  5.0  X  10®  psi  =  0.75  x  10®  psi 

E  =  1.5  X  10®  psi  =  0.25 

zz  ^  Lz 

The  interlaminar  epoxy  resin  phase  has  the  properties: 

E  =  0.5  X  10®  psi  V  =  0.35 

m  ^  m 

G  =  0.185  X  10®  psi 
m  ^ 
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Interlaminar  Crack  Growth  During  Fatigue 

Microscopic  observation  in  the  fatigue  experiments  revealed  that 
delamination  was  initiated  after  a  few  load  cycles.  It  propagated  in  a 
stable  manner  and  was  clearly  a  matrix-dominated,  progressive  failure 
mechanism  occuring  in  a  resin-rich  region  between  plies.  This  suggested 
an  appropriate  mechanical  model  used  later  in  the  analytical  study.  Mea¬ 
surements  were  made  on  the  change  of  the  subcritical  crack  length  during 
cyclic  loading.  Typical  fatigue  crack  growth  data  in  unidirectional  and 
(±45°)^^  glass/epoxy  composites  are  presented  in  Figs.  5  and  6.  Test 
results  shown  in  the  figures  for  several  maximum  cyclic  stress  levels 
indicate  that  crack  lengths  change  drastically  with  the  number  of  load 
cycles  and  that  crack  growth  rates  are  very  distinct  between  the  unidi¬ 
rectional  and  angle-ply  laminates.  The  growth  rates  of  the  delamination 
d£^/dN  were  determined  by  a  computational  procedure  based  on  the  first 
derivatives  of  crack  growth  curves  fitted  by  experimental  data.  The  crack 
extension  rates  were  found  to  increase  continuously  during  the  course 
of  fatigue,  and  rapid  crack  propagation  occurred  when  the  ultimate  delamination 
resistance  of  the  composite  was  approached.  Intralaminar  failure,  which 
may  lead  to  a  complete  fracture  of  the  composites,  were  often  observed  in 
(45°/-45°//  45°/-45°/-45'’/45°/-45°/45°)  and  quasi-isotroplc  (45°/-45°// 
0°/90°/90°/0°/-45°/45°)  glass/epoxy  prior  to  large-scale  delamination. 

The  number  of  loading  cycles  to  fracture  versus  applied  cyclic  stress 
for  the  two  composite  systems  is  shown  in  Fig.  7.  The  ultimate  static 
strengths  of  these  laminates  were  14.7  and  44.55  Ksi,  respectively. 


12 


Interlaminar  Crack-Tip  Response 

Bie  driving  force  of  the  interlaminar  crack  growth  depends  primarily  on 
local  crack-tip  stresses.  Thus,  it  is  necessary  to  have  accurate  infor¬ 
mation  on  the  crack-tip  behavior.  Since  the  geometry  and  material  proper¬ 
ties  are  very  complex  in  the  problem,  the  response  of  the  composite  to  the 
damage  may  be  best  illustrated  by  isostress  contours.  These  contours  were 
plotted  from  the  analytical  solutions  obtained  by  the  aforementioned 
hybrid  finite  element  analysis .  Figure  8  gives  the  distribution  of  the 

in-plane  stress  a  /a  around  the  crack  with  5, ,  =  0.8  inch  in  a  (0°/0°// 

XX  «>  d  - 

0°/0°/0°/0°/0°/0°)  glass/epoxy.  The  computer  plot  provides  a  graphic 
description  of  the  solutions  and  an  overall  view  of  the  damage  behavior. 

It  reveals  several  salient  features  of  the  problem:  (i)  the  intensified 
stress  field  caused  by  interlaminar  damage  is  very  localized  within  the  ma¬ 
trix-rich  interlaminar  region  ahead  of  the  crack  tip,  (ii)  a  large  stress 
gradient  exists  through  the  laminate  thickness  with  high  concentration  in 
adjacent  plies,  and  (iii)  a  local  compressive  in-plane  stress  component 
is  developed  along  the  lower  crack  flank.  Figures  9  and  10  show  the  contour 
plots  of  interlaminar  shear  and  transverse  normal  stress  distributions  in  the 
near  field.  Both  stress  components  reach  extremely  high  levels  within  the 
thin  interlaminar  region  when  the  crack  tip  is  approached.  The  interlaminar 
stresses  extend  continuously  through  the  interface  into  the  adjacent  plies. 

A  more  accurate  description  of  the  near  field  stresses  may  be  expressed 
by  a  log  versus  log  r  representation,  where  r  is  the  distance  ahead  of 
the  crack  tip.  The  numerical  results  near  the  crack  tip  give  straight  lines 
with  a  slope  of  -i,  indicating  a  classical  fracture  mechanics  solution,  i.e., 
a  l//r  singularity;  details  of  the  results  may  be  found  in  [l5  1  •  The  opening 
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and  shearing  mode  stress  components  are  observed  to  have  the  same  order  of 
magnitude  at  the  interlaminar  crack  tip  when  a  uniaxial  nominal  stress  is 
applied.  The  singular  stress  domain  is  rather  small  and  is  embedded  in 
the  resin-rich  interlayer  region. 

The  redistribution  of  laminate  stresses  due  to  the  propagating  inter¬ 
laminar  crack  will  have  a  significant  effect  on  the  deformation  and  fracture 
of  the  composite.  For  example,  the  compressive  stress  in  the  lower  crack 
flank  may  Introduce  severe  local  buckling  of  the  delaminated  plies  in  the 
vicinity  of  the  surface  notch.  This  phenomenon  was  experimentally  observed 
in  [16,17].  The  highly  nonuniform  stress  distributions  in  the  transverse 
direction  shown  in  the  isostress  contour  plots  suggest  that  the  classical 
laminate  theory  may  not  be  applicable  to  a  delaminated  composite.  The  stress 
concentration  may  also  Introduce  a  transition  from  interlaminar  to  intra¬ 
laminar  fracture  mode  when  different  ply  orientations  are  used  [18].  Further¬ 
more,  the  presence  of  the  interlaminar  crack  may  become  detrimental  to  the 
serviceability  of  the  material  when  used  in  an  adverse  environment  such  as 
mositure  and  corrosive  chemicals,  because  not  only  does  the  crack  expose  the 
interior  of  the  composite  to  environmental  attack,  but  also  the  stress  concen¬ 
tration  in  the  near  field  may  enhance  the  rate -depen dent  degradation  process, 
which  is  evident  by  the  results  in  [19]. 

Fatigue  Crack  Growth  Rate  Characterization 

The  stable  extension  of  an  interlaminar  crack  in  glass/epoxy  composites 
reveal  several  features  of  the  progressive  interlaminar  fatigue  damage ,  which 
can  be  used  to  probe  the  nature  of  the  delamination  growth  behavior.  Since 
the  interlaminar  crack  growth  is  observed  as  a  matrix -dominated  failureprocess 
and  since  it  is  strongly  affected  by  the  local  crack-tip  response,  it  may  be 
appropriate  to  postulate,  as  in  the  cases  of  metals  and  polymers[20,21] , 
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that  AKj  and  are  critical  controlling  factors  in  the  fatigue  damage. 

The  ranges  of  the  cyclic  stress  intensity  factors  are,  in  turn,  related  to 
the  crack  length,  the  cyclic  stress,  the  geometry  and  laminate  properties. 

In  the  current  fatigue  test,  stress  intensity  factors  change  continuously 
as  the  load  and  the  interlaminar  crack  length  increase,  and  the  corre¬ 
sponding  crack  growth  rate  varies  accordingly.  Fora  given  crack  length 
and  laminate  configuration  ,  exact  stress  intensity  factors  may  be  computed 
by  using  the  present  analytical  method.  Therefore,  AK^.  and  can  be 

precisely  determined  from  Eqs.  9  and  10.  It  may  be  reasonable,  at  this 
point,  to  postulate  further  that  AK^.  and  independently  affect  the 

interlaminar  crack  growth  rate.  Thus,  dil^/dN,  obtained  from  the  experiments, 
may  be  related  to  both  the  AK^  and  the  simultaneously  as  shown  in 

Figs.  11  and  12  for  the  glass/epoxy  composites  with  two  distinct  ply  orienta¬ 
tions.  Strong  correlation  between  dil^/dN  and  AKs  is  observed  from  the 
results,  and  a  power  law  relationship  of  the  following  form  may  be  established. 


log 

d 

dN 

log(AK^)  +  log(AKj^)  +  C2 

“l  “2  “3 


(11) 


where  a.s  are  the  directional  cosines  of  the  lines  determined  by  a  least- 
1 

square  curve-fitting  method,  and  and  are  constants  associated  with 
the  threshold  values  of  the  mixed-mode  stress  intensity  factor  ranges . 
This  functional  relationship  in  Eq.  11  suggests  that 


or 

where  the  exponents 
crack  geometry  and 


d£  X 


(12a) 


d£  A 

'V  (AKj.^) 


(12b) 


A^  and  A^  are  related  to  the  laminate  configuration,  the 
the  material  properties  of  composites.  The  results  in 
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Figs.  11  and  12  give  respective  and  X2  the  values  of  7.32  and 
7.42,  for  0°  and  7.52;  and  7.74  for  ±  45°.  The  strong  correla¬ 
tion  between  the  interlaminar  crack  propagation  rate  and  the  ranges  of 
cyclic  stress  intensity  factors  are  very  informative.  Further  study  is 
needed  to  identify  the  nature  of  the  interlaminar  crack  under  more  complex 
fatigue  loading  and  to  construct  the  crack  growth  rate  surface  of  mixed- 
mode  failure  so  as  to  characterize  fully  the  interlaminar  fatigue  behavior,, 
and  the  relative  dominance  of  each  individual  fracture  mode. 

Limitations  of  Current  Study 

It  should  be  noted  here  that  several  geometric  and  material  compli¬ 
cations  in  composites  have  not  been  included  in  the  current  analytical 
modeling,  such  as  the  nonuniformity  of  the  resin-rich  interlaminar  region 
thickness,  the  heterogeneity  of  the  fiber  and  matrix  phases,  and  other 
types  of  commonly  encountered  defects  (e.g.  voids,  inclusions,  broken 
plies,  etc.).  These  complications  are  statistical  in  nature  and  almost 
always  exist  in  the  material.  Their  presence  may  have  significant  effects 
on  the  interlaminar  crack  growth  behavior,  but  has  not  been  fully  explored 
yet.  The  anti-plane  mode  deformation  and  fracture  (mode  III),  which  may 
have  an  important  contribution  to  the  progressive  fatigue  failure  of  the 
materials,  especially  in  the  angle— plied  composite  laminates,  have  not  been 
considered  in  the  present  study.  A  recent  study  on  edge  delimination  [10] 
reveals  that  the  mode  III  component  of  the  crack— tip  field  has  a  very  high 
value,  which  may  be  primarily  responsible  for  the  interlaminar  crack  growth 
and  failure  of  the  composites,  including  the  l  45°  case.  Furthermore,  resi¬ 
dual  stresses  due  to  ply  differential  thermal  contraction  during  the  curing 

process  may  strongly  affect  the  crack  growth  behavior.  This  effect  has  not 
been  included  in  the  present  investigation. 
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SUMMARY  AND  CONCLUSION 

A  combined  experimental  and  analytical  investigation  of  interlaminar 
crack  growth  in  fiber-reinforced  composites  under  cyclic  loading  has  been 
presented.  The  unidirectional  ply  configuration  was  studied  first  to 
elucidate  the  fundamental  failure  behavior  and,  later,  served  as  a  reference 
for  studying  more  complex  angle-ply  composite  laminates .  In  the  experimental 
phase,  microscopic  observation  of  failure  mechanisms  revealed  that  delami¬ 
nation  is  clearly  a  matrix-dominated,  progressive  failure  occurring  in  the 
resin-rich  interlaminar  region.  Interlaminar  crack  growth  rates  under 
tension-tension  fatigue  were  obtained,  and  test  results  were  further  used 
in  conjunction  with  analytical  characterization.  The  analytical  study  of 
the  problem  was  based  on  an  advanced  hybrid-stress  finite  element  method, 
formulated  by  using  Muskhellshvili's  complex  stress  functions  through  a 
modified  complementary  energy  principle.  Numerical  solutions  clarify  the  basic 
nature  of  the  cyclic  interlaminar  failure.  Implicaticns  of  the  results  in 
the  applicability  of  classical  laminate  theory  to  the  delamination  problem 
due  to  the  localized  interlaminar  crack-tip  singularity  and  the  stress  con¬ 
centration  in  the  adjacent  plies  were  discussed.  Fatigue  crack  growth  rates 
were  found  to  be  related  directly  to  the  ranges  of  the  mixed-mode 
cyclic  stress  intensity  factors  by  a  power  law  relationship.  The  study 
elucidates  many  important  features  of  the  delamination  problem,  provides 
insight  into  the  interlaminar  fatigue  behavior,  and  should  be  of  practical 
Importance  in  the  testing  and  design  of  composites  materials. 
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MICROGRAPH  OF  GRAPHITE  FIBER  REINFORCED 
COMPOSITE  SHOWING  INDIVIDUAL  PLIES  AND 
RESIN -RICH  INTERLAYERS  BETWEEN  PLIES. 


a.)  BEFORE  LOADING 
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b.)  DURING  LOADING 

FIG.  2  FATIGUE  TEST  SPECIMEN  FOR  INTERLAMINAR  CRACK 
GROWTH  STUDY  (2L'=  6.0  in.  ;  W=I.O  in.  ) 
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STRESS  vs.  TIME  FOR  FATIGUE  TEST 


FIG.  4  INTERLAMINAR  CRACK  GEOMETRY,  COORDINATE  SYSTEM  AND 

HYBRID-STRESS  F.E.M.  MESH  CONFIGURATION  IN  COMPOSITES 
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FIG.  5  INTERLAMINAR  CRACK  GROWTH  DURING  FATIGUE 
I N  (QW  W°/070°/07070°)  GLASS/  EPOXY 
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FIG.6  INTERLAMINAR  CRACK  GROWTH  DURING  FATIGUE  IN 
(45%45W5%457-45°/45%45745°)  GLASS /EPOXY 
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FIG.  7  FATIGUE  STRESS  vs  NUMBER  OF  CYCLES  TO  FAILURE  FOR  GLASS/EPOXY 
COMPOSITES  WITH  DIFFERENT  LAMINATE  CONFIGURATIONS 
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FIG. 8  IN-PLANE  LONGITUDINAL  STRESS  CONTOUR  akx/C&,  NEAR 

INTERLAMINAR  CRACK  TIP  IN  (0°/D°//070y0*n)707D°) 
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